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Abstract. The GOV star HD 166435 has been observed by the fiber-fed spectrograph ELODIE as one of the 
targets in the large extra-solar planet survey that we are conducting at the Observatory of Haute-Provence. We 
detected coherent, low-amplitude, radial-velocity variations with a period of 3.7987 days, suggesting a possible 
close-in planetary companion. Subsequently, we initiated a series of high-precision photometric observations to 
search for possible planetary transits and an additional series of CaH H and K observations to measure the level 
of surface magnetic activity and to look for possible rotational modulation. Surprisingly, we found the star to be 
photometrically variable and magnetically active. A detailed study of the phase stability of the radial-velocity 
signal revealed that the radial- velocity variability remains coherent only for durations of about 30 days. Analysis 
of the time variation of the spectroscopic line profiles using line bisectors revealed a correlation between radial 
velocity and line-bisector orientation. All of these observations, along with a one-quarter cycle phase shift between 
the photometric and the radial- velocity variationss, are well explained by the presence of dark photospheric spots 
on HD 166435. We conclude that the radial- velocity variations are not due to gravitational interaction with an 
orbiting planet but, instead, originate from line-profile changes stemming from star spots on the surface of the 
star. The quasi-coherence of the radial-velocity signal over more than two years, which allowed a fair fit with a 
binary model, makes the stability of this star unusual among other active stars. It suggests a stable magnetic field 
orientation where spots are always generated at about the same location on the surface of the star. 
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1. Introduction 

A large extra-solar planet survey has been under- 
way since 1994 at the Observatory of Haute-Provence 
with the high-precision, fiber-fed echelle spectrograph 
ELODIE ( Baranne et al. 1996| ) mounted on the 193cm 
telescope. The search for extra-solar planets is carried 
out by seeking changes in the radial velocity of each 
star produced by gravitational interaction with orbit- 
ing planets. In 1995 this program lead to the first de- 
tection of an cxtrasolar planet orbiting a Sun-like star 
( Mayor fc Qucloz 1995| ). Our survey contains 324 G dwarf 
stars brighter than V — 7.65, including HD 166435 
( iQucloz ct al. 1998b| ). 

The star HD 166435 is a GO dwarf with (B-V)= 0.633. 
From its distance of 25 pc and its apparent magnitude 



my = 6.85, we find an absolute magnitude My = 4.8. 
If the metallicity is roughly solar, HD 166435's loca- 
tion in the H-R diagram suggests a main-sequence age of 
a few billion years. Interestingly, the star was detected 
in the ROSAT all-sky survey as an extreme ultravio- 
let (EUV) source, but no other evidence of activity was 
found by MuUiss & Bopp (1994) in their analysis of the 
Hq and Call (8542A) lines. [Burleigh et al. 1997| found a 
hint of emission in the CIV (1549 A) line but detected no 
other obvious spectral emission features. When we se- 
lected HD 166435 for ELODIE observations, we had no 
real indication that HD 166435 might be young and ac- 
tive. A possible coincidence with another EUV source had 
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been mentioned (Burleigh et al. 1997) as a possible expla- 
nation for the EUV ROSAT detection. 

After only a few radial-velocity measurements of 
HD 166435, it became apparent that the star exhib- 
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ited low-amplitude, radial-velocity variations with a pe- 
riod of about 4 days. Further measurements in the 
following months confirmed our initial finding, lead- 
ing us to believe that HD 166435 may have a plane- 
tary companion in a close orbit. Subsequently, we be- 
gan a campaign of high-precision photometric observa- 
tions with an automatic photoelectric telescope (APT) 
at Fairborn Observatory in Arizona to search for plane- 
tary transits and made tentative plans to announce the 
new planet at the Protostars and Planets IV confer- 



ence (Mannings, Boss, & Russell (2000)). To our surprise 



however, we found the star to be photometrically variable 
with the same period as the radial-velocity variations. In 
this article we describe these data, along with additional 
Call H and K measurements acquired at Mount Wilson 
Observatory, and discuss the most likely explanation for 
the observed variations in all three data sets. 

2. ELODIE spectra 

We obtained 70 spectra with ELODIE from 1997 May 
to 1999 September. These spectra have typical signal- 
to-noise ratios of 70-150 (at 5000A) per resolution el- 
ement and a resolution (A/AA) of about 42,000. The 
data reduction was carried out on-line during the observa- 
tions by the ELODIE automatic reduction software (see 
Baranne et al. 1996 for details). All the spectra were ac- 
quired in the high-precision mode, which provides a si- 
multaneous thorium reference spectrum. The zero-point 
wavelength calibration has an intrinsic instrumental pre- 
cision of 10ms~^. 

To examine certain key spectral features, we co-added 
the 70 individual spectra to obtain a composite spectrum 
with a very high signal-to-noise ratio of about 1000. Since 
the individual spectra were sampled at slightly different 
velocities relative to the star due to the earth's orbital 
motion, the composite spectrum also has a better quality 
than the individual spectra. Two key spectral features in 
the composite spectrum are shown in Fig. 1. The Call H 
spectral region shows an emission reversal in the core of 
the absorption line, suggesting strong chromospheric ac- 
tivity. The broader Ca II photospheric wings compared to 
the Sun are the consequence of the higher effective tem- 
perature of HD 166435. No strong lithium absorption line 
is seen. The comparison with the Sun's spectrum suggests 
that, with our resolution, the sensitivity to the ^Li doublet 
is limited to medium or strong features. Moreover, in the 
case of HD 166435, the Lithium doublet is blended with a 
nearby feature due to rotational broadening. However, it 
appears that the ''Li feature in HD 166435 is no stronger 
than the one in the solar spectrum. We estimate an upper 
limit of 10 mA for the ^Li doublet at 6707.8A . The curve 



of growth from Soderblom et al. (1993a) then implies an 
upper limit for the lithium abundance of N(Li)< 1.7. 

We computed the radial velocity of each individual 
spectrum with the ELODIE automatic reduction soft- 
ware, which uses a cross-correlation technique with binary 
mask templates. In addition to the stellar radial velocity. 
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Fig. 1. Selected regions of the composite spectra of 
HD 166435 (solid line). For com parison, a solar integrated 
spectrum ( Kurucz et al. 1984 ) with the same resolution 
has been superimposed (hatched line). Upper panel 
Call H line. An emission feature is clearly visible in the 
core of the line. Lower panel Region of the Li^ line at 
6708A . The solar spectrum is slightly offset upwards for 
display purpose. No strong Li^ feature is detected in the 
HD 166435 spectrum. Notice in both figures the larger line 
broadening of HD 166435 compared to the Sun 



the cross-correlation function (CCF) provides additional 
information about features in the averaged spectral lines 
(see Queloz 1994 for a review). The width of the CCF 
yields the v sini of the star, while the equivalent width of 
the CCF can be used as a metallicity estimate if the tem- 



perature of the star is known approximately (Mayor 1980 
Bcnz fc Mayor 1981 1 ; [Queloz 199^ ). The mean equivalent 



width of the CCF for the 70 individual spectra is similar 
to the value measured for the Pleiades stars of the same 
temperature (B— Vq = 0.63). Therefore, a solar metallicity 
can be assumed for HD 1664 35. With the same t echnique 
and calibration described in [Queloz et al. 199^ we der- 
vive a mean vsuii= 7.6 ± 0.5 km s^"'^ . A similar value is 
measured by CORAVEL {vsmi= 7.7 ± 0.7kms"i). 

2.1. Radial-velocity data 

Analysis of our entire two-year radial- velocity data set re- 
vealed a periodicity of 3.7987(±.0004) days. The radial 
velocities are shown in Fig. 2 where they are phase folded 
on that period. We fit these data with a binary model and 
derived a slight orbital eccentricy e = 0.2 and an ampli- 
tude K — 83 ms^^, corresponding to a minimum mass for 
a planetary companion of 0.6 Mj. The rms of the residuals 
to the fit is 28 ms""'^. A large value compared to the typi- 
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Fig. 3. Variation of the rms (0-C) of a sine-curve fit to 
the radial-velocity data versus the time duration of the 
observations. 



Fig. 2. Phase diagram of the radial-velocity data with a 
period of 3.7987 day. A binary model is superimposed on 
the data (solid line). See text for details. 

cal ELODIE precision of 10 ms^^. If the binary model is 
the correct interpretation of the observed radial- velocity 
variations, then additional radial- velocity scatter intrinsic 
to the atmosphere of the star must be present. 

In order to test the binary hypothesis, we first investi- 
gated the coherence time of the radial-velocity variations. 
Since a short-period planet in a (presumed) circular or- 
bit induces sinusoidal reflex motion in the star, we fit sine 
curves with a period of 3.7987 days to subsets of the radial- 
velocity observations of increasing duration and computed 
the rms of the residuals in each case. The results of this 
test are shown in Fig. 3. The rms of the sinusoidal fits in- 
crease with an increase of the duration of the data sets. For 
durations loss than 30 days, the residuals are in agreement 
with the expected rms by the photon and instrumental 
noise. For durations longer than 30 days, phase disruptions 
likely occur and limit the phase coherence and therefore 
increase the rms of the fit. 

To quantify the phase disruptions in the radial- velocity 
data, we searched for phase changes among 21-day blocks 
of data spread over the two-year duration of our observa- 
tions. From Fig. 3 above, we sec that, for data sets of three 
weeks duration, the radial- velocity variations are coherent 
enough to be described as a sinusoid to within the preci- 
sion of the measurements. The top panel of Fig. 4 shows 
phase shifts of up to ± 0.1 cycle among the 21-day blocks 
of data. For comparison, we applied the same procedure to 
the 51 Peg radial- velocity data, where the circular orbital 
motion of this confirmed planet in close orbit produces si- 
nusoidal reflex motion. As expected, no significant phase 
changes can be seen in the 51 Peg data (Fig. 4, bottom 
panel) . Interestingly, the random phase changes in the ra- 
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Fig. 4. Phase shifts in the radial velocities derived by fit- 
ting sine curves to subsets of the data with a maximum 
duration of 21 days, (a) HD 166435 data, (b) 51 Peg data. 
Note that there is an arbitrary time offset between the two 
data sets for display purposes. The two hatched lines in 
each panel indicate the zblcr error on the global fit. 

dial velocities from HD 166435, while signifcantly larger 

than in 51 Peg, are not so large as to prevent their reason- 
ably good description by a sinuosidal model of the whole 
two-year data set. 

The detection of significant phase shifts in the 
HD 166435 radial velocities explains the high rms 
(28ms~^) of the residuals to the planet model fit in Fig. 2. 
We also examined the amplitude K of the fits and their 
zero levels Vq for similar changes with the same technique. 
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No significant change of the K amphtude was detected for 
either star. The Vq values for HD 166435 exhibited small 
variations that are probably related to line-profile varia- 
tions (see below). 

2.2. Line profile and bisector analysis 

The radial velocity of a star is defined to be the velocity 
of the center of mass of the star along our line of sight. 
However, the observational determination of a star's radial 
velocity is accomplished by measuring the Doppler shift 
of spectral lines produced in the stellar atmosphere. If 
changes occur in the star's spectral-line profiles, the mea- 
sured Doppler shifts may not correspond precisely to the 
velocity of the star's center of mass. For pulsating stars 
and stars with very active atmospheres, this is indeed the 
case. 

To interpret observed radial- velocity variations as true 
changes in the velocity of the star, we must show that the 
observed variations do not stem from changes in the stellar 
atmosphere. One of the best ways to do this is to look for 
subtle changes in the line profiles. Radial- velocity changes 
in a star accompanied by spectral-line profiles that are 
constant in time would be a direct evidence for real ve- 
locity changes in the star's center of mass. The detection 
of line-profile changes would point instead to atmospheric 
phenomenon. 

The computation of a stellar radial velocity from a par- 
ticular spectrum with the CCF technique involves the cor- 
relation of thousands of spectral lines. Thus, in order for 
line-profile changes to affect the radial-velocity measure- 
ments, the profile changes in a particular spectrum must 
affect all of the spectral lines in a similar way. Otherwise, 
no change would be induced in the cross-correlation func- 
tion. 

Our CCF technique employs a binary template 



( iQucloz 1994| ) and is a very efficient way to look at global 
phenomena in spectra. The CCF function represents a 
mean spectral-line profile of all of the lines selected by 
the template. To study the mean spectral-line profile of 
HD 166435, we designed a specific template that selects 
only the weak and non-saturated spectral lines. We used 
the line-bisector technique to analyze the CCF profiles in 
the same way that we would analyze an individual spec- 
tral line. In Fig. 5 we show HD 166435's mean CCF and 
its bisector. The bisector displays the classic "C" shape 
observed in other stars lying on the same (cool) side 



of the granulation boundary (Gray 1989). Unlike similar 
studies ( [Toner fc Gray 1988 ), our data have an absolute 
wavelength calibration with an accuracy of approximately 
10ms~^ (the precision of the instrument). Therefore, we 
don't need to shift the bisector arbitrarily to a mean value 
for comparison. Moreover, if a bisector change is detected, 
we know whether it is the top or the bottom of the bisector 
that is moving. 

We selected two regions of the bisector at the top and 
bottom of the profile {Vt and Vf,) (see Fig. 5 for illustra- 
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Fig. 5. Right The mean CCF function of HD 166435's 
spectra constructed with a template selecting only the 
weak and non-saturated lines. This profile represents the 
mean spectral-line profile of the lines selected by the tem- 
plate. Left The bisector of the CCF. Vq is an arbitrary 
offset. Note the definition of the boundaries for the com- 
putation of (Vt and Vf,). 

tion) to measure possible bisector orientation changes with 
high accuracy. The difference between Vt and Vf, is equiva- 
lent to the " bisector velocity spa n" widely used in similar 
studies (see Toner fc Gray 198^ ). This is actually a mea- 
surement of the inverse of the mean slope of the bisector. 

Any correlation between radial-velocity changes and 
line-bisector orientation would cast serious doubts on the 
reflex-motion interpretation of the radial-velocity varia- 
tions. If the radial- velocity variations are due to center-of- 
mass velocity changes of the star, we would see the bisec- 
tors oscillating to the left and right of the mean bisector 
with no change in shape or orientation. In order to investi- 
gate this, we first looked at the bisectors of each individual 
CCF. In Fig. 6 the bisectors of the CCFs for two sets of 
spectra selected at opposite phases of the radial-velocity 
cycle {(p — 0.0 ± 0.1 and (j) — 0.5 ± 0.1) are displayed. 
Comparison with the mean bisector shows that the bi- 
sectors are twisting around the mean bisector with the 
same periodicity as the radial-velocity signal. The top of 
the CCF profile moves very little, while the bottom shifts 
back and forth with an amplitude 50% larger than the K 
amplitude of the radial-velocity signal. The only way to 
understand this twisting motion of the bisectors is through 
a change in the spectral-line profiles synchronized with the 
radial- velocity cycle. The amplitude of this effect is enough 
to produce the observed radial-velocity signal. 

To investigate further the exact relationship between 
the orientation of the bisectors (bisector span) and the ra- 
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Fig. 6. Individual bisectors for two sets of spectra selected 
at opposite phases of the radial-velocity cycle, a Spectra 
measured at ip — 0.0 ± 0.1. b Spectra measured at cj) — 
0.5 ± 0.1. The hatched line illustrates the mean bisector 
computed by averaging all spectra. 

dial velocity of the CCF, we plot these two parameters 
in Fig. 7. A direct relationship between the two quan- 
tities is clearly visible. A linear solution (Vt — Vb) = 
-0.88(±0.04)K. can be fit with a 27 ms^^ rms. Thus, the 
results of our line-profile analysis of HD 166435's CCFs 
provide strong evidence that the radial- velocity variations 
originate in the stellar atmosphere and not from reflex 
motion of the whole star. 



3. Photometry 

The Hipparcos catalogue lists 144 photometric measure- 
ments of HD 166435 acquired between 1990 January and 
1993 March and identifies the star as an "unsolved" vari- 



able with a range of 0.05 mag (Ferryman et al. 1997). Our 
periodogram analysis of the Hipparcos photometry con- 
firms the lack of periodicity, particularly near the radial- 
velocity period. This is what led us initially to suspect 
that our observed radial-velocity variations might be due 
to the presence of a planet in an orbit similar to that of 



the companion of 51 Pegasi (Mayor & Queloz 1995). 

Thus, shortly after the detection of the periodicity 
in the radial-velocity signal, we began photometric ob- 
servations of HD 166435 with the T8 0.80 m automatic 
photoelectric telescope (APT) at Fairborn Observatory in 
Arizona to search for possible transits of the companion 
across the disk of the star. We acquired 326 Stromgren b 
and y observations with a two-channel precision photome- 
ter on the APT between 1998 June and 2000 June. The 
observations were reduced differentially with respect to 
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Fig. 7. Radial velocity of each CCF versus the bisector 
span (Vt — Vb) of the CCF profile. The dotted fine is the 
best linear fit to the data. 



the comparison star 99 Her (HR6775, HD 165908, F7V), 
corrected for atmospheric extinction with nightly extinc- 
tion coefficients, and transformed to the Stromgren sys- 
tem. External precision of a single observation with the 
0.80 m APT averages 0.0011 mag. Further details of the 
automatic-telescope operations and data-reduction proce- 
dures can be found in Henry (1999)| . The individual pho- 
tometric observations are available on the Tennessee State 
University Automated Astronomy Group web site[|. 

Periodogram analysis of the entire set of 326 ob- 
servations taken together reveals a photometric period 
of 3.7995 ± 0.0005 d. Thus, the photometric and radial- 
velocity periods agree within their respective uncertain- 
ties. 

For our photometric analysis we adopt, as refer- 
ence, the radial-velocity period of 3.7987-day and To — 
2450996.5, the radial-velocity maximum of the best-fit 
sine-curve model. The photometric observations are di- 
vided into five groups as shown in Table 1, where we give 
the results of least-squares, sine-curve fits on the radial- 
velocity period to the five Stromgren y data sets. We 
also list the periods derived from the five individual data 
sets. The photometric amplitudes vary from 0.050 mag to 
0.035 mag; the mean magnitudes of the five data sets have 
a range of about 0.023 mag. Although the five light curves 
approximate sinusoids, the rms values are, nevertheless, 
somewhat higher than the 0.0011 mag precision of typical 
observations. This is primarily the effect of slight cycle-to- 
cycle changes in the light curves within each data set; the 
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Table 1. Summary of the Stromgren y Photometric Results 



Data 


Date Range 


Nobs 


Period 


Mean Brightness 


Full Amplitude 


Phase of 


rms 


Set 


(JD - 2400000) 




(days) 


(mag) 


(mag) 


Minimum 


(mag) 


1 


50986-50996 


219 


3.732(6) 


1.849(1) 


0.050(1) 


0.20(1) 


0.003 


2 


51086-51136 


16 


3.776(7) 


1.865(1) 


0.052(4) 


0.12(1) 


0.004 


3 


51224-51360 


42 


3.807(10) 


1.857(1) 


0.025(4) 


0.07(3) 


0.008 


4 


51434-51501 


25 


3.763(34) 


1.842(2) 


0.023(4) 


0.23(3) 


0.007 


5 


51591-51710 


24 


3.821(7) 


1.842(2) 


0.032(6) 


0.29(3) 


0.008 



rms values increase for the later data sets since the later 
data sets are longer and contain more cycles. 

Data set 1, which covers only 10 days or 2.7 cycles, has 
the largest number of observations since at this time we 
were making repeated measurements each night to search 
for transits. These data are plotted in the left and right 
third panels from the top of Fig. 8. Based on the Hipparcos 
photometry mentioned above, we did not expect to find 
coherent light variability on this period. Instead, we found 
a smoothly varying, nearly sinusoidal light curve with a 
period matching the radial- velocity period and a minimum 
near phase 0.25, where we had hoped to find transits. 

In Fig. 8 we see that the time of radial-velocity max- 
imum, i.e., maximum red shift at phase 0.0, corresponds 
closely to the time when the light curve crosses its mean 
level at Ay = 1.85 mag in the downward direction. The 
time of photometric minimum corresponds closely to the 
time when the radial-velocity curve crosses its mean level 
(zero velocity shift with respect to the mean velocity) in 
the downward direction. The time of radial-velocity mini- 
mum corresponds closely to the time when the light curve 
crosses its mean level in an upward direction. Finally, 
the light-curve maximum corresponds closely to the time 
when the radial-velocity curve crosses its mean level in 
an upward direction. Hence, these two nearly sinusoidal 
curves are one-quarter cycle out of phase with each other. 

Our Stromgren photometry and A{b — y) colors from 
data set 1 are plotted in the bottom two panels of 
Fig. 8, also with phases computed from the radial- velocity 
ephemeris. The color curve is exactly in phase with 
the light curve and has a peak-to-peak amplitude of 
0.009 mag. The star is redder when it is fainter. 

4. Ca II H and K Spectrophotometry 

From 1998 June to 1998 August, 141 Ca II H and K mea- 
surements were made on 20 nights with the 100-inch tele- 
scope at Mount Wilson Obs ervatory as part of the HK 
Project ( Baliunas et al. 1998 ). In that program, measure- 
ments of the Ca II H and K lines of several thousand stars 
are made as a proxy for surface magnetism. 

The observed quantity, S, is the fiux measured in two 
0.1 nm pass bands centered on the H and K lines normal- 
ized by two 2.0 nm-wide sections of photospheric fiux cen- 
tered at 390.1 nm and 400. Inm. A nightly calibration fac- 
tor is determined from measurements of a standard lamp 
and standard stars ( Baliunas et al. 1995 ). The night-to- 



night rms precision of the lamp is on the order of 0.5%, 
while the standard stars have an average standard devia- 
tion of ~ 1.5%, which limits the lowest amplitude of vari- 
ability that can be detected to approximately 1%. 

Our periodogram analysis of the 141 Ca II observa- 
tions over a period range of two to ten days gives a strong 
periodicity at 3.85 ± 0.01 d, essentially identical to the 
radial-velocity and photometric periods. Since the time 
span of the Ca II observations is only 62 days and the 
period uncertainty is larger than for the radial velocities, 
we adopted the radial-velocity period and Tq = 2450996.5 
for the analysis. 

A 30-day subset of the S measurements, during which 
we had simultaneous photometry and radial- velocity data, 
is plotted in the second from the top panels of Fig. 8. 
A sinusoidal fit to these S values gives an amplitude of 
0.023 and a mean value of 0.46. Interestingly, the phase 
is shifted 1/8 of a cycle compared to the light and color 
curves (with 15% statistical error). The variation of S is 
therefore 1/8 cycle out of phase with the radial- velocity 
variation. The star has its higher level of Ca II H and K 
emission 1/8 of a cycle after it is faintest and 1/8 of a 
cycle before it reaches its most red-shifted radial velocity. 
Similar phase differences between photometric and pho- 
tospher ic data are also d etected in the GO Hyades star 
VB31 (iRadick et al. 1987|). 



5. Discussion and modeling 

Our observations of HD 166435 indicate that the star is 
active and thus young. The relatively fast rotation also 
indicates a young star. If we assume that the X-ray fiux 



detected by ROSAT ( [Burleigh et al. 1997| ) indeed origi- 
nates from the star, the 0.5 c/s PSPC c ounts in the 0.1- 
2.4keV range ( Randich fc Schmitt 1995 ) corresponds to 
log(Lx) = 29.5 or log(ix/iboi) = -4.1. About the 
same X-ray luminosity is measured in other single stars 
of similar temperature in the 220 Myr cluster NGC 6475 
( James fc Jeffries 1997 ). Thus, we estimage the age of 
HD 166435 to be about 200 Myr. 

The typical Li content of GO stars in the Pleiades 
is about log(A^Li) w 3.0. At the age of the Hyades 
(400 Myr), the lithium is depleted to about 2.6. The so- 
lar value is 1.0. The weak lithium content of HD 166435 
noted above (log(A^Li) < 1.7) conflicts with its young 
age deduced from its rotation rate and activity level. 
We have no explanation for this. However, other young 




Fig. 8. Left Simultaneous observations of (from the top) radial velocity, S index, delta y magnitude, and delta [b — y) 
color, of HD 166435 over a time span of 30 days. A best-fit sine-curve with a period fixed at 3.798 days is shown. A 
vertical line is drawn at an epoch of maximum radial velocity to help visualize the phase offsets between data sets. 
Right Same data but phase folded with P = 3.798 d and To = 2450996.5. 



stars with low lithium content have been detected in 
Preasepe (KW 392) ([Soderblom ct al 1993bD and Hyades 
(vB9, vB143) ( frhorburn et al. 19931 ). 

The brightness of HD 166435, its spectroscopic line 
profiles, and its chromospheric emission all vary with the 
3.8-day periodicity found in the radial velocities. As seen 
in Fig. 8 above, the brightness variations exhibit a roughly 
one-quater phase shift with respect to the radial veloci- 
ties. All of this strongly suggests that the radial-velocity 
variations are due to the observed stellar activity. For a 
magnetically active star with dark spots being carried into 
and out of view by rotation, we would expect to see such 
a phase shift between the brightness changes and radial- 
velocity variations. When the spot lies on the central 
meridian of the star, it lies partly on the approaching and 
partly on the receeding halves of the stellar disk. The spot 
also causes the maximum depression in the brightness of 
the star at this time. Hence the radial- velocity effect of the 
spot at photometric minimum should be null. As can be 
seen in Fig. 8 above, the radial-velocity shift is very close 



to zero at the minimum of the light curve. Radial- velocity 
minimum (the time of maximum blue shift) should oc- 
cur as the dark spot approaches the receeding limb of the 
star and hence is moving out of view. This results in more 
light from the approaching half of the stellar disk reach- 
ing our telescopes than from the spotted, receeding half 
of the disk. This is in agreement with Fig. 8, where we see 
that radial- velocity minimum corresponds to a time when 
the star is rapidly brightening. Similar reasoning predicts 
that radial- velocity maximum should occur when the star 
is rapidly dimming, as also observed in Fig. 8. The delta 
{b — y) color curve in the bottom panel of the figure shows 
that the star is redder when it is fainter, further showing 
that the spots causing the brightness variation must be 
cooler (darker) than the surrounding photosphere. 

The small phase shift between the brightness of 
HD 166435 and its chromospheric emission level suggests 
that the photocenter of the dark spot distribution is some- 
what offset from the photocenter of the bright plage area. 
This phenomena is observed on the Sun to a lesser degree 
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where leader spots tend to be larger and to live longer 
than follower spots, resulting in a displacement up to 5" 
between the photocenters of the spot and the associated 
plages ( Dobson-Hockey fc Radick 1986| ). The larger dis- 
placement in HD 166435 suggests that the size of active 
regions is much larger than those seen on the Sun, in agree- 
ment with the observed amplitude of the light curve. 

The growth and decay of individual spots on time 
scales longer than a stellar rotation will affect the phase 
and amplitude of the radial-velocity signal (as well as 
the brightness and chromospheric emission levels) on long 
time scales. Our detection of coherent radial-velocity vari- 
ations without noticeable phase shifts over intervals of 
up to 21 days suggests that the spots are stable for at 
least this long. This provides an estimate of the evolu- 
tionary time scale for individual spots of roughly a month. 
However, we observed a quasi-coherent signal for the en- 
tire two-year span of our radial- velocity measurements. 
This suggests that new spot formation occurs at about the 
same longitude on the surface of the star during this time. 
Among the numerous active stars measured in our extra- 
solar planet survey, the stability of the radial- velocity sig- 
nal of HD 166435 is unique. None of the other active stars 
in our sample display the level of phase coherence observed 
in HD 166435 that could suggest the possibihty of plane- 
tary reflex motion. Therefore, the spot generation mecha- 
nism of HD 166435 may be the result of unusual stability 
in the geometry of its magnetic field. This phenomenon is 



not unique amongst active stars. Toner & Gray 1988 ob- 



served a stable phase lasting for at least 3 years in the 
variation of the line bissector of the young star x Boo A. 
Stable active regions at the same longitude are observed 



as well in RS CVs binaries ( Jetsu 1996). Moreover, on the 



Sun, during flares, hot spots have been se en for 30 years 
at the same longitude ( Bai 1988 ; Bai 1990| ) 

We conclude from these qualitative arguments that a 
rotating spot model can successfully explain the variability 
in our spectroscopic and photometric data sets. Moreover, 
the mean level of the Ca index, < S >= 0.46, correspond- 
ing to R'ffj^ = —4.26, is in a good agreement with the 
activity level expected for a star with rotation period of 
3.8 days (|Noyes et al. 1984| ). 

The combination of our rotation period and 
V sin i measurements allows us to estimate that the in- 
clination of the star's rotation axis to our line of sight is 
approximately i = 30'^. Notice that in our case the in- 
crease of the macroturbulence with the star rotation is 
negligible. It does not change our estimate of the i an- 
gle. The smooth (i.e., without plateaus) variations in our 
radial velocities, photometry, and chromospheric emission 
suggest that the largest spotted region is always visible to 
the observer and does not completely disappear when on 
the far side of the star. Therefore, it is likely that the spot 
region is located within 30 degrees of the stellar pole and 
is visible, at least in part, throughout the rotation cycle. 
In this case, projection effects would be responsible for 
most of the rotational modulation of the star's brightness. 



Thus, we have shown that surface magnetic activity 
on the star HD 166435 can mimic the kind of radial- 
velocity variations observed in stars with true planetary 
reflex motions. However, further analysis of the spectr- 
copic line profiles, photometry, and chromospheric emis- 
sion clearly demonstrates that stellar activity is the ori- 
gin of the radial-velocity variations. Our observations of 
this peculiar object actually strengthen the interpretation 
of other low-mass companions to solar-type stars where a 
low-amplitude, radial-velocity variation has been detected 
but where photometric, c hromospheric, and line-profile 



variations are absent (e.g., Henry et al. 2000 ) 
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